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The two Higgs doublet model has a rich vacuum structure, including the possibility of existence 
of two Standard Model-like minima at tree-level. It is therefore possible that the universe's vacuum 
is metastable, and a deeper minimum exists. We present the analytical conditions one must demand 
of the potential's parameters to prevent that possibility, and analyse what the current LHC data 
tells us about the eventual existence of that second minimum. 

I. INTRODUCTION 

I The two Higgs doublet model (2HDM) 1J is one of the simplest extensions of the Standard Model (SM) of 

particle physics, in which the number of scalar doublets is doubled compared to the SM. This simple addition 



makes for a richer scalar spectrum, which includes two CP-even scalars, the lightest h and the heaviest H , 

T± 



a pseudoscalar, A, and a charged scalar, H ± . The 2HDM boasts very interesting phenomenology, including 
possible spontaneous CP violation, tree-level flavour changing neutral currents mediated by scalars and dark 
matter candidates. For a recent review, see [2J. The recent discovery of the Higgs boson [5J 0] enables us to 
further constrain the 2HDM parameter space, and ascertain whether the model survives comparison with data. 
In fact, the 2HDM does a very good job describing the LHC results [SHE]- 

In the SM there is only one possible type of vacuum. In the 2HDM, on the other hand, charge breaking vacua 
may occur, and in those situations the photon would acquire a mass. Also, as was already mentioned, we also 
have possible minima in which the CP symmetry is spontaneously broken, alongside with the electroweak gauge 
symmetry. But a remarkable property of the 2HDM potential is that these different types of minima cannot 
coexist |10[ lllj . We call a vacuum which breaks electroweak symmetry but preserves the electromagnetic and 
CP symmetries a "normal" minimum. And whenever a normal minimum exists, any possible charge or CP 
breaking stationary points are necessarily saddle points, and lie above the normal minimum. The stability of 
the normal minimum against charge or CP breaking is thus guaranteed. But there is still another difference 
in the 2HDM vacuum structure regarding the SM: for certain choices of parameters, the 2HDM potential can 
have two normal minima simultaneously [121414] . The two Higgs doublets, $1 and $2, would have vacuum 
expectation values (vevs) v\ and V2 such that vf + v| = 246 GeV 2 in one of the minima, and all elementary 
particles would have the known masses - this would be the minimum where the universe is currently at. In the 
second minimum, however - and this minimum could be deeper or higher than ours - the fields would have vevs 
{ v i > w 2}j with v'1 + v'2 7^ 246 GeV 2 - and the elementary particles might have masses much smaller, or larger, 
than what is observed. 

It is therefore possible that "our" vacuum is not the state of lowest energy of the theory. Thus the universe 
would be in a metastable state, with the possibility of tunneling to the deeper vacuum. We call this situation 
the "panic vacuum" . In [TS] we presented the conditions that the parameters of the potential need to obey so 
one can avoid the presence of a panic vacuum in the softly broken Peccei-Quinn [TB] version of the 2HDM. We 
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also showed that the LHC data already excludes most of the parameter space where panic vacua might occur 
in this model. In this talk I will discuss the existence of double neutral minima in the 2HDM potential with a 
Z 2 discrete symmetry. That is the most used version of the 2HDM, and panic vacua do occur in it, for plenty 
of the model's parameter space. The bounds one has to impose on the potential to avoid such minima were 
deduced in [5] and will be reviewed here, as well as what the current LHC data tells us about the existence of 
the second minimum. 



II. THE VACUUM STRUCTURE OF THE 2HDM 



The most general 2HDM potential has 14 real parameters, which can be reduced to 11 using the reparametriza- 
tion invariances of the model. That model, however, leas to tree-level flavour changing neutral currents mediated 
by scalars. In order to avoid them - since they are incredibly constrained by experimental measurements - one 
usually imposes a discrete Z 2 symmetry on the lagrangian, such that <!>i — > <&i and $2 — > — $2- This was 
first proposed by Glashow, Weinberg and Paschos [El [18], and the resulting scalar potential has only seven 
independent real parameters. Its parameter space is however extremely constrained, which is why one many 
times adds a softly breaking term m\ 2 . The softly broken potential is therefore written as 



V = ml^ + ml^ 



4 ($1 



$9 + h.c. 



~Ai|$i| 4 + * A 2 |$ 2 | 4 + A 3 |$i| 2 |$ 2 | 2 + A 4 |$ t 1 $ 2 | 2 + ~A 5 



h.c. 



(1) 



where we have also imposed CP symmetry on the scalar potential and all the parameters shown are real. 
As we have already mentioned, the 2HDM can have charge breaking vacua, wherein the doublets acquire 
vevs ($i)cb = (0, ci) T /y/2 and ($ 2 )cb = ( c 2 , c^) T /v2- If, simultaneously, there is a "normal" solution 
(vi , v 2 )/y/2 of the minimization equations, the depth of the potential at the charge breaking stationary point, 
Vcb, is related to the depth of the potential at the normal stationary point, Vn, by [TtH IPT] 



Vcb — Vn — 



Av 2 



[(V1C3 



\2 1 2 21 
V 2 Ci) 



(2) 



A' 



where rn 2 H± is the square of the charged scalar mass at the normal stationary point - meaning, if that is a 
minimum we will have Vcb — Vn > 0; further, it was also shown that the charge breaking stationary point 
is, in this case, necessarily a saddle point. Ergo, the normal minimum is the global one, and no tunneling 
to a deeper minimum can occur. Likewise, if there is a CP breaking stationary point, with vevs such as 
(&i)cp = (0 , vi) T '/v2 and ($ 2 )cp = (0 , v 2 e lS ) T /\/2, the difference of depths of the potential at this stationary 
point and a normal one is given by 



Vcp — Vn — 



\Av 2 



) [(v 2 v\ cos 9 — viv 2 ) 2 + v^y 2 sin 2 9] 

) N 



(3) 



where m 2 A is the pseudoscalar mass at the normal stationary point. Again, if the normal stationary point is a 
minimum, it will be the deepest and thus stable against tunneling. In short, the existence of a normal minimum 
guarantees that the global minimum of the potential is also normal. However, the 2HDM, due to the soft 
breaking term, can have two normal minima [T2TU4j . Further, it may be shown that if the potential has a depth 
equal to Vn for the minimum with vevs {v± , v 2 }, and a depth Vn> for the minimum with vevs {v[ , v' 2 }, we 
have 



V N > - V N 



N 

V1V2 
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In this expression, both the charged Higgs mass m 2 H± and the sum of the squared vevs v 2 is evaluated at each 
minima. It is not therefore clear which is the deepest minimum. But in [TH1 HH HH] the conditions for the 
existence of two neutral minima in the 2HDM were established. They were cast into a simpler form in [T5] for 
the Peccei-Quinn model, and in [9] for the Z 2 one. In [9] we also presented a thorough deduction of the bounds 
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to avoid panic vacua, and their generalisation for all 2HDM CP-conserving potentials. Therein we were able to 
deduce simple necessary conditions for the existence of two normal minima for any such potentials, which we 
will not reproduce here for briefness. Fortunately, the study of panic vacua does not require that one analyses 
whether or not the potential satisfies those conditions. In fact, all we need do is compute the following quantity 
D, which we call a discriminant: 

D = m^m 2 ^ — fc 2 ™ 2 ^) (tan /3 — fc) with k— \ '/-r— • (5) 

V A 2 

As usual, tan j3 = v^/vx, written with the vevs of "our" vacuum. The existence of a panic vacuum is thus 
summarised in the following theorem: 

Our vacuum is the global minimum of the potential if and only if D > 0. (6) 

It is therefore extremely simple to ascertain whether the 2HDM minimum is, or is not, the global minimum 
of the model: all one has to do is to compute the value of D above, having reconstructed the parameters of the 
model from experiments (this may prove to be difficult, but it is, in principle, quite achievable). Up to date, 
there is no evidence of the existence of any scalars beyond that which is predicted by the SM. Nonetheless, as 
we will now see, the current LHC data already allow us to exclude most of the parameter space where panic 
vacua might occur. 



III. PANIC VACUA EXCLUSION USING LHC DATA 



The 2HDM scalar potential we are considering has 8 independent parameters. That is a vast parameter space, 
but it can be constrained in many different ways. In the first place it needs to have one minimum with vevs 
which originate the correct masses for the W and Z bosons, i.e. the doublets must need have vevs v\ and V2 
such that vf + v 2 = 246 GeV 2 . Further, since the discovery of the Higgs boson, we must ensure that the lightest 
CP-even scalar has a mass of 125 GeV. On the other hand, the potential has to be bounded from below, which 
imposes the following constraints on its quartic couplings: 

Ai > , A 2 > , 

A3 > -AAiAa , A 3 + A 4 - |A 5 | > -V^iAa ■ (7) 

Then, just as in the SM, we need to ensure that the model obeys perturbative unitarity, which again constrains 
the potential's quartic couplings [20l [21]. And of course, the model also must comply with the precision 
electroweak constraints, via the bounds on the S, T and U parameters [2"2"H2"5"| . However, all of these bounds 
apply exclusively to the scalar sector of the theory, but one must be aware that there are plenty of restrictions 
on extended scalar sectors from their interactions with fermions, namely from B-physics experiments. The 
Z2 symmetry imposed on the scalar potential is extended to the fermion sector in a variety of ways, we will 
consider only the two most studied: in model Type I the fermion fields transform under the global Z 2 in such 
a way that only $ 2 couples to the fermions; and in model Type II 4> 2 couples only to up-type quarks, and $1 
to down-type quarks and charged leptons. Both of these models possess very different phenomcnologies, and 
B-physics bounds have been taken into account in our simulations (26] [27] . With all of these constraints in 
place, we performed a vast scan of the 2HDM parameter space, for both models considered, taking rrih = 125 
GeV, 125 < m H < 900 GeV, 90 < m A ,m H ± < 900 GeV, -tt/2 < a < tt/2, 1 < tan/3 < 40 and \m\ 2 \ < 900 
GeV 2 . At this point we are ready to compute the observables which will be compared to the LHC data. 

In order to reduce the large uncertainties present in the calculation of hadronic cross sections at the LHC, 
we consider the ratios between the observed rates of the Higgs boson decaying into certain particles and their 
expected SM values. If we assume that what is being observed is described by the 2HDM, that ratio is defined, 
for a given final decay state / of the lightest CP-even Higgs boson h, as 

_ a 2HDM {pp -> h) BR 2HDM (h -» /) 
/_ a SM {pp~>h)BR SM (h^ f) ' ( ) 

Thus both production cross sections a and branching ratios (BR) of the Higgs boson are included, and we have 
considered all possible Higgs production mechanisms at the LHC. 

In Fig. [T] we present, for both models under consideration, the rates of the light Higgs h into two Z bosons 
versus the rate of h into two photons. In green (light grey) are all the points obtained in our scan, with the 
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Model I Model II 




FIG. 1: Rzz as a function of i? 77 for models I and II. Green (light) points represent all points generated, red (dark) ones 
correspond to the existence of a panic vacuum (D < 0). The solid (dashed) lines correspond to la (2a) experimental 
bounds on the rates. 

above constraints, both theoretical and experimental. In red (black) we show the points for which a panic 
vacuum occurs - meaning, points for which the value of D, calculated from eq. ([5|, is negative. It is important 
to consider that the density of points generated is so large that there are many green points scattered in the 
middle of the red ones - meaning, the areas marked red are not necessarily excluded (a consequence, of course, of 
the fact that we are dealing with an 8-dimensional parameter space, and these figures are only 2-dimensional) . 
Fig. ^clearly shows that there are regions, in the plane R 11 -Rzz, which are completely free from panic vacua. 
The solid and dashed lines shown in the plots correspond to conservative la and 2a intervals on the combined 
values for Rzz and i? 77 , Rzz = 0.93 ±0.28, i? 77 = 1.66 ±0.33, which we took from ref. [35], based on the LHC 
data before the Moriond conference [2"5] . Notice that after the recent Moriond update on the LHC results [30 
these numbers may have changed substantially, but there isn't yet an official combination of the ATLAS and 
CMS results. However, the plots we show in this communication have the advantage of being easily adapted 
for changing LHC experimental bounds, by simply drawing over them different black lines. 

The remarkable thing is how much the current LHC data already can tell us about the nature of the 2HDM 
vacuum, even if no extra scalars have been found. In fact, as can be seen from Fig. [I] the panic points are 
distant from even the 2a bands, which include some non-panic region as well, for model Type I. That does not 
occur in model Type II, in which some of the panic region is inside the 2a region. But notice that there are 
many choices of parameter space values still allowed by the current data for Model II which do not lead to panic 
vacua. Thus, at least in these two variables, both types of models are capable of describing the current data. 
Nonetheless, that data does not exclude the possibility, in model type II, of our vacuum being metastable. 

It has been possible to measure at the LHC - with considerable uncertainty - the production of Higgs bosons 
via different processes, namely gluon-gluon fusion and vector boson fusion (VBF). Analysing these processes 
separately gives us information about the coupling of the Higgs to both fermion (the gluon-gluon process) and 
gauge bosons (the VBF one). We shall use the results of the ATLAS experiment [28], which appear as la and 
2a ellipses in the R^-R^ F plane. Our results appear in Fig. [2J for both Type I and Type II models. We 
observe that the experimental exclusion of points with panic vacua is not as thorough as that which happened 
with the previous observables. In model Type I it is not possible to exclude, at 2a, the existence of panic vacua. 
However, the panic vacua points which now seem allowed have been excluded in Fig. [I] For model Type II, 
even the la bands include panic vacua solutions. We observe, nonetheless, that the ellipses contain plenty of 
green/light grey points as well - which means that there are many allowed choices of parameters for which panic 
vacua do not occur. We even see that in these variables the Type II model agrees with the data at the la level, 
something which the Type I model cannot achieve. 

The current results for a Higgs decaying to rf are compatible with the expected SM value. In fact, ATLAS 
measured R TT = 0.7 ± 0.7 and CMS R TT — 0.72 ± 0.52. We can see how the panic vacua points are distributed 
in the {R TT , i? 77 } in Fig. [3] The rr data (we represent the less restrictive bounds, those of ATLAS), tells us 
that panic vacua are 2a disfavored in model Type I, and that model barely agrees, at 2a, with the LHC results 
in Rjj, agreeing at la in R TT )- In model Type II, there are many panic vacua solutions not excluded by the 
data at la; but for much of Model II's parameter space, we have agreement with the experimental results at 
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FIG. 2: i? 77 , with Higgs production via gluon-gluon fusion versus vector boson fusion production, for models I and II. 
Same colour codes as Fig. [T] 
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FIG. 3: R TT as a function of R 11 for models Type I and Type II. Same colour codes as Fig. [T] 



the ltr level, with or without panic vacua. 



IV. CONCLUSIONS 



The rich vacuum structure of the 2HDM includes the possibility of two neutral minima being able to coexist 
at different depths. Thus there is the possibility that the vacuum we are currently living in is metastable, which 
we called the panic vacuum. It is possible to find an extremely simple criterium, the discriminant D of eq. ^ 
being positive, which, if obeyed, guarantees that no metastability occurs. We emphasise that this situation 
is quite different from the SM one: there a metastable (or even unstable) vacuum may develop but only due 
to radiative corrections. In fact, the importance of radiative corrections to the results found here cannot be 
overstated, and remains an open question. We also performed an estimate of the lifetime of the false vacuum 
in the panic situation, and verified that it is, for the vast majority of the points scanned, much inferior to the 
age of the universe - as such, these panic vacua are indeed to be excluded, since "our" vacuum would have 
decayed long ago. Nevertheless, we see that the current LHC results already tell us a lot about the stability 
of the vacuum in the 2HDM. For instance, a measurement of Rzz and i? 77 very close to 1, with sufficient 
precision, would exclude the possibility of panic vacua. We also saw that the values of the parameters of the 
potential which produce panic vacua do not correspond to uninteresting regions of the model - rather, they 
predict observables which are not absurd and indeed may fall into the current experimental bounds. This, by 
itself, indicates the need to take seriously this possibility of vacuum instability in the 2HDM. 
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